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Abstract. In this paper we construct an "Interior Penalty" Discontinuous Galerkin method 
to approximate the minimizer of a variational problem related to the p(x)— Laplacian. The 
function p : SI — > [pi-.pz] is log Holder continuous and 1 < p\ < p-z < oo. We prove that the 
1 minimizers of the discrete functional converge to the solution. We also make some numerical 

experiments in dimension one to compare this method with the Conforming Galerkin Method, 
in the case where pi is close to one. This example is motivated by its applications to image 
processing. 



< 

J3 1 1- Introduction 

+-> • 

In this paper we study a discontinuous Galerking method to approximate the minimizer of a 
non homogeneous functional. This functional is related to the so-called p(x)— Laplacian operator, 
i.e. 

: (1.1) A p{x) u = div(\Vu(x)\P^- 2 Vu). 

\Q . This operator extends the classical Laplacian yp(x) = 2) and the so-called p— Laplacian 

(p(x) = p with 1 < p < oo) and it has been recently used in image processing and in the 
^! ■ modeling of electrorheological fluids, see [3, 8, 23] . 

Ov 

In an image processing problem, the aim is to recover the real image I from an observed 
image £ of the form £ = / + r], where rj is a noise. 

Approaches to image denoising have been developed along three main lines: wavelet methods, 
stochastic methods and variational methods, see references in [3]. One variational approach that 
has attracted a great deal of attention is the total variation method of L. Rudin, S. Osher and 
E. Fetami [22]. The variational problem is 



Minimize the functional |Du|(0) over all the functions in BV{Q) n L 2 (J7) such that 

and /|«-{|^ = „> 
Jn Jn Jn 

for some a > 0. 

The conditions on the space come from the assumption that £ is a function that represents a 
white noise with mean zero and variance a. Moreover, the authors prove that this problem is 
equivalent to minimizing 



\Du\(Q) + - [ \u-i\ 2 dx 
2 Jn 
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for some nonnegative Lagrange multiplier A = A(cr, £). This model works when the image is 
piecewise constant, but in some cases can cause a staircasing effect. See for instance [7]. 
An older approach consists in minimizing, 



This method solves the staircasing effect, but it has the problem that it does not preserve edges. 

In [8], the authors introduce a model that involves the p(x)— Laplacian, for some function 
p : fi — >■ [pi , 2] , with pi > 1. This function encodes the information on the regiones where the 
gradient is sufficiently large (at edges) and where the gradient is close to zero (in homogeneous 
regions). In this manner, the model avoids the staircasing effect still preserving the edges. 

Recently, in [3] the authors propose a variant of the method of Chen, Levine and Rao [8]. 
More precisely, they consider the functional 



with p : £1 — > [1, 2] a function such that p(x) = Pm(\VGs * £\(x)), where G§(x) is an approxima- 
tion of the identity, M » 1 and Pm is a function that satisfies Pm(0) = 2 and Pm(x) = 1 for 
all \x\ > M. Observe that, since p(x) = 1 for some values of x, the authors have to rewrite the 
functional in a form that allow for computation of weak derivatives. 

Motivated by the above mentioned applications, we study a numerical method to approximate 
minimizers of a functional related to the p(x)— Laplacian. 
We work in the following setting: 

Let SI be a bounded Lipschitz domain. For functions p,s,t the following conditions will be 
assumed when necessary, 

(HI) p : O — > [pi,p2] (1 < Pi < P2 < oo) is log-Holder continuous. That is, there exists a 
constant C\ og such that 






(H2) s G L°°(n), with 1 < s(x) < p*(x) - e for some e > 0; 
(H3) t G C°(dn) with 1 < t(x) < p*(x). 



Here, p* and p* are the Sobolev critical exponents for these spaces, i.e. 



(1.2) 




and 




Given p,q,r satisfying (H1),(H2) and (H3) respectively and £ E L q ('\Q) 
the functional 



we want to minimize 




over all v € A, where 



A = {v G W 



v — ud G W, 



(«)}, 
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ud € W 1 '^'' (fl) and <90 = Td UT^r. For the definitions of the variable exponent Sobolev spaces 
W l ^\Q) and Wp^ ( ) (0), see Section 2. 

Let us observe that, considering the applications we have in mind, it is relevant to study 
the minimization problem in the case where p approaches the value 1 in some regions. We can 
see, by making some numerical experiments, that the minimizers have large derivative in these 
regions. For this reason, the Conforming Finite Element Method is not appropriate since, its 
use would imply the need of fine meshes in order to obtain good approximations, see Section 8. 

We consider the so-called Discontinuous Galerkin Methods. These methods are relatively new 
from the theoretical point of view. In [1], we can find a unification of all methods of this type. 
In all the examples of that paper, the authors take as model a linear differential equation. 

At this point we want to mention that, in [3] and [8] the authors find an approximation of 
the solutions by using an explicit finite difference scheme for the associated parabolic problem. 

Our aim is to study, in the future, the minimization problem in the case where p approaches 
the value 1 in some regions (where there is no weak formulation). For this reason, we think 
that the best way to find approximations is by finding a good discretization of the minimization 
problem. We take a discretization similar to the one in [6] where the authors study a functional 
that includes the case p = constant. 

Our discrete functional is the following 

h(v h ) = jf (\Vv h + R h {v h )\ p{x) + \v h - dx + jf \v h - un^h 1 -^ dS 

+ / |Ml p(l) h 1 ~ p(a;) d5+ f \v h \ r( > x US, 

where h is the local mesh size, h is the global mesh size, is the union of the interior edges of 
the elements, [Ju/j]] is the jump of the function between two edges and Vt> h denotes the element- 
wise gradient of v^, see Section 3 for a precise definition. Observe that the boundary condition 
is weakly imposed by the second term of the functional. Lastly, is the lifting operator defined 
in Section 5.1, which represents the contributions of the jumps to the distributional gradient. 

In the case p = constant, the boundedness of this operator is proved by using an inf-sup 
condition. In our case, for technical reasons, we use a different approach, which consists of 
finding a local characterization of the operator. 

With this setting the discrete problem is to find a minimizer Uh of 1^ over the space S k (J~h) 
of all the functions that are polynomials of degree at most k in each element, with k > 1, see 
Section 3. 

In this work we show that the sequence Uh converges to the minimizer u oi I over the space 
A. We want to remark here that we are assuming that, for each Vh € S k {Th), all the terms of 
Ih{vh) can be exactly computed. 

In fact, we prove the following 

Theorem 1.1. Let be a polyhedral domain. Let p(x) , q(x) , and r(x) be functions satisfying 
(HI), (H2) and (H3) respectively and let u D G W 2 ' P2 (n). For each h 6 (0,1], let u h G S k (%) 
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be the minimizer of Ih- If u is the minimizer of I then 



(1.3) Uh — > u strongly in L s ' '(f2) Vs satisfying (H2), 

(1.4) Uh — > u strongly in L^'^dQ) Vi satisfying (H3), 

(1.5) JfcK) -> J(u), 

(1.6) #(u h ) ->■ strongly in L p ^(n), 

(1.7) / |u h -u D | p ^h 1 - p W ( i5+ / IWI^h^^^^O, 

(1.8) Vu h ^ Vu strongly in L p{ '\tt). 



Lastly, we want to mention the places where we need the regularity hypotheses on the func- 
tion p. First, in order to prove Theorem 1.1 we need to use the continuity of the embedding 
W lj *>(fi) ^ L p *( )(fi), the continuity of the Trace operator W 1 ^^!) ^ L^(dn) and the 
Poincare inequality. As we can see in Theorem 2.9, Theorem 2.10 and Theorem 2.8 for these 
results we need p to be log- Holder and r € C°(<9S1). 

We also use strongly that p is log- Holder in Proposition 2.11. This result says that if n is an 
element with diameter h K and p+ and pi are respectively the maximum and minimum of p over 

k then h P K ~ P+ is bounded independent of h K . This property is crucial in the proof of several 
results along the paper. 

On the other hand, in order to prove the convergence of the sequence Uh we need a technical 
hypothesis on the boundary condition up. In fact, Lemma 6.4 only covers the case where 
u D e W 2 > P2 (Q). 

Outline of the paper. In Section 2 we state several properties of the Variable Exponent 
Sobolev Spaces. 

In Section 3 we give some definitions and properties related to the mesh and to the Broken 
Sobolev Spaces. 

In Section 4 we study the Reconstruction operator and we prove some error estimates that 
are crucial for the rest of the paper (Corollary 4.5). 

In Section 5 we prove the boundedness of the Lifting operator (Theorem 5.3). 

In Section 6 we prove the Broken Poincare inequality (Theorem 6.1), the coercivity of the 
functional (Theorem 6.2) and finally we give the proof of Theorem 1.1. 

In Section 7 we study the convergence of the Conforming Finite Element Method. 

In Section 8 we give a Id example and compare both conforming and non-conforming schemes. 

2. Preliminaries: The spaces L p ^'\Q) and W 1,p ^(Q) 
We now introduce the spaces LfW(O) and W^Xty and state some of their properties. 

Let p: f2 —> \pi,P2] be a measurable bounded function, called a variable exponent on £1 where 
Pi := essinf p(x) and p 2 := esssupp{x) with 1 < p\ < p 2 < oo. 

We define the variable exponent Lebesgue space L P ^'\Q) to consist of all measurable functions 
u : — > M. for which the modular 

P (.)W := / Hx)f x Ux 
Jn 
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is finite. We define the Luxemburg norm on this space by 

IMIz,p(-)(n) = IMU-) := [nf { k > 0: e P (-)( u / k ) ^ !}■ 

This norm makes L p ('\fl) a Banach space. 

The following properties can be obtained directly from the definition of the norm, 
Proposition 2.1. Ifu,u n S L p ('\$l), ||it||p(.) = A, then 

(1) A < 1 (= 1, > 1) if only if [ \u{x)\ p{x) dx < 1 (= 1, > 1); 

Jo 

(2) IfX > 1, tfien A Pl < / jn(x)| p(x) < A P2 ; 

(3) //A < 1, then X P2 < [ \u(x)\ p( - x) dx < A Pl ; 

Jn 

(4) / \u n {x)\ p{x) dx^O if only if \\u n \\ p{ .) -> 0; 



(5) ||l|| p (.) < max||n|"i , | ; 



are open sets then there exists a constant C > depending 

m 

|MLi>(-)(Q) < I Ml £*>(•) fO:V 



8=1 

on m suc/i that 



i=l 



□ 



Proof. See Theorem 1.3 and Theorem 1.4 in [18]. 

For the proofs of the following three theorems we refer the reader to [21]. 
Theorem 2.2. Let q(x) < p(x), then L p ('\Q) L q( -'\il) continuously. 

Theorem 2.3. Let p,q,r : £1 — >• [1, oo) and e > be such that p(x) < r(x) < q(x) — e for all 
x £ Q.. Then, there exists a positive constant C such that for every u E L P ^(J7) n the 
inequality 



\u\\ L r(. )(Q) <C\\uf T , )(n A\ U 



holds, where fj> > and v > are define as 



esssupn 



essinf n 



v = < 



essinfc 



p(x 




— r(x) 


r(x] 


qix) 


— p(x) 


p{x) 


q(x) 


— r{x) 


r(x) 


q(x) 


— p{x) 


<?(X 


r(x] 


— p{x) 


r(a^ 


q(x) 


— p{x) 


q(x) 


r(x) 


-p{x) 


r(x) 


q(x) 


-p{x) 



if 
if 
if 
if 



L«<» (fi) 

u \\lp(-)(Q) > 1) 

u \\Lp(-)(n) ^ 1' 

u llz<r(-)(n) > 1) 

u llz,9(-)(n) - 



Theorem 2.4. Letp'(x) such that, l/p(x) + l/p'(x) = 1. Then L p ' ^(Sl) is the dual of L p{ -\£1). 
Moreover, if pi > 1, L p ('-*(f]) and are reflexive. 
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Now we give some well known inequalities, 
Proposition 2.5. For any x fixed we have the following inequalities 

\V ~ e {x) < C(\ V \^- 2 V - \H\ P ^- 2 0(V - ifp(x) > 2, 

\V ~ Z\ 2 (\V\ + \Z\Y {X) ~ 2 < C(\r,\^- 2 r, - \e {x) ~ 2 Z)(v ~ ifp(x) < 2, 



|^| P (x) < 2^)- 1 (|r ? - £\P(*) + if p [x) > 1. 

These inequalities say that the function A(x,q) = \q\ p ( x ' 2 q is strictly monotone. 
Proposition 2.6. Let F n ,F G L P ^(Q). 
(1) V 



then 

(2) // 
then 



F n -^ F weakly in LpW(O) 

/ |F| p ^ds < liminf / |F n | p ^dx. 
n ^°° Jn 

F n F strongly in 

[ \F n \P^dx^ ( \F\P^dx. 
Jn Jn 



(3) // 

(2.9) F n ^ F weakly in L P ^(Q) and [ \F n \ p{x) dx ^ [ \F\ p{x) dx 

Jn Jn 



then 



F n ^ F strongly in L P ^(Q). 



Proof. For the proof of (1) and (3) see Theorem 3.9 and Lemma 2.4.17 in [13]. Finally (2) follows 
by Proposition 2.3 in [15]. □ 

Let W^itt) denote the space of measurable functions u such that, u and the distributional 
derivative Vu are in LPO(fi). The norm 

IMIiyi>p(-)(n) : = IMIp(0 + lll Vu lllp(-) 

makes W l,p<y '\d) a Banach space. 

We define the space Wq ,p{) {Q) as the closure of Cg°(n) in WMOffi). Then we have the 
following version of Poincare inequality, see Theorem 3.10 in [21]. 

Lemma 2.7. If p : O — > [1, +oo) is continuous in Q, there exists a constant C such that for 
every u E Wo' pi '\n), 

\\u\\lp(-)(0) < C\\Vu\\ LP (.)( n y 

We also have the following version of the Poincare inequality, see Lemma 2.1 in [20], 

Theorem 2.8. Let C M. n be a Lipschitz domain and p, q : CI — > [l,+oo) with p < q < p* . 
Then, 

\\u - («)n|| L ?(-)(n) < C\\Vu\\ LPi . ){Q) 
for all u 6 W 1,P ^'\VL), where (u)q = pr f n udx. 
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In order to have better properties of these spaces, we need more hypotheses on the regularity 
of p. For expample, it was proved in [11], Theorem 3.7, that if one assumes that dSl is Lipschitz 
and p is log-Holder continuous then C°°(ft) is dense in W lj "^(Q), see also [10, 14, 16, 21, 24]. 
The local log-H16der condition was first used in the variable exponent context in [26]. 

We now state two Sobolev embedding Theorems (for the proofs see [12] and Corollary 2.4 in 
[17], respectively). 

Theorem 2.9. Let Q be a Lipschitz domain. Letp : — > [1, oo) be log-Holder continuous. Then 
the embedding W l ' p ('(Fl) '—t L p *('\fl) is continuous. 

Theorem 2.10. Let Q be an open bounded domain with Lipschitz boundary. Suppose that 
p € C°(0) with pi > 1. If r € C°(dQ) satisfies the condition 1 < r(x) < p*(x) for all x € dCl, 
then there is a compact boundary trace embedding W 1,p ^(n) L r ('\d£l). 

Let Td C dfl, and p be log-Holder. We define the space \ty as the closure of the space 

{(p E C°°(n) : ip = on F D } in W 1 ^)^). 

The following proposition is crucial in order to prove the main result of this paper. 

Proposition 2.11. Let p : — > [l,oo) be log-Holder continuous and bounded. Let a > 0, 
DcO and h = diam(D) then: 

(1) There exists a constants C independent of h such that 
(2.10) h a(p(x)- P (y)) < c VX)y G rj ; 

(2) If A>h a then A^ x ) < CA^ for all x,y £~D such that p(x) < p(y). 

Proof. Let x,y £ D. If p(x) > p(y) or h > 1 the result follows since fl is bounded. If p{x) < p{y) 
and h < 1, using that p is log-Holder, we have 

p(y) ~ P(x) < — 7 — ^— ; — \- < 



log I e + L~^i J lo § ( e + X 

Then, we get (2.10). 

By (2.10) and as A > we have that for all x,y € D such that p(x) < p(y), 

/ a \ p{%)-p(y) 

j^>(x) _ [ _ I i l a(p{x)-p(y)) < (Jj\p(y) 

h° 



□ 



3. The mesh 7^ and properties of W 1,p (')(7/i) 

In this section we give some definitions and properties related to the mesh and to the Broken 
Sobolev Space. 

Hypothesis 3.1. Let Q be a polygonal Lipschitz domain and (Th)he(o,i] be a family of partitions 
of £1 into polyhedral elements. We assume that there exist a finite number of reference polyhedral 
«i, k r such that for all k £ 7~h there exists an invertible affine map F K such that, k = F K (ki). 
We assume that each k E Th is closed and that diam(K) < h for all k E Th- 
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Now we give some notation, 

£ h ■= { K n k' : dxm H (K n k') = N - 1} U {k D dfi : dim#(K n 50) = iV - 1}, 
Imi := (J{e G £ h : dim H (e D dfl) < N - 1}. 

Mh is the set of nodes of Th- For every z G Ah and e € 4 we define, 

T z := |J{Ase 7^: z€k}, T k := J{T Z : z G k}, T e := \J{T K : e G «}, 
/i K := diam(K), /i 2 := diam(T 2 ) /i e := diam(e), 

:= ess inf p(x) p+ := esssupp(x), pi := essinf p(x) and p+ := esssupp(x). 

iGk irG/t xC = e xGe 

We assume that the mesh satisfies the following hypotheses, 

Hypothesis 3.2. The family of partitions (Th)he(o,i] satisfies the Hypothesis 3.1 and 

(a) There exist positive constants C\ and Ci, independent of h, such that for each element 

K G 7h 

CiC < kl < c 2 h%. 

(b) There exists a constant C\ > swc/i £aa£ /or a// /i G (0, 1] and /or all face e £ £h there 
exists a point x e G e and a radius p e > C\diam{e) such that B Pe (x e ) PiA e Ce, where A e 
is the affine hyperplane spanned by e. Moreover, there are positive constants such that 

ch K <h e < Ch K , ch K i <h e < Ch K > 

where e = kCi k' . 

We use the notation ~ to compare quantities which differ only up to positive constants that 
do not depend on the local or global mesh size or on any function which appears in the estimate. 

Remark 3.3. By the regularity assumption of the mesh we have the following, 

B{ze^:zGK}~l, §{k £%■ k C T z } ~ 1, 

tf{«' G %■■ «' C T K } ~ 1, j}{e G £ h : e C T z } ~ 1 and f{e6^:eC T K } ~ 1. 

Remark 3.4. As a consequence, we have that diam(T K ) ~ h K and for each z G k and e C 5k, 
^2 ~ h K and /i e ~ h K . See the discussion on Section 4.2 in [6]. 

Remark 3.5. By Proposition 2.11, we also have that for each edge e C 8k, ~ /ie f° r an y 
£ k. We will replace p1,p e _ by p_ and p+,p+ by p+ when no confusion can arise. 

Now, we introduce the finite element spaces associated with Th- We define the variable broken 
Sobolev space as 

W lj, U(%) ■= {u G L\n): u\ K G W 1A -\k) for all k G T h }, 

and the subspaces 

U k {T h ) := {u G C(Q) : «|„ G P fc for all k G 7~4, 

^(Tfe) := {u G L J (n) : u\ K G P fc for all k G 
where P k is the space of polynomials functions of degree at most k > 1. 
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We also define, for any k G Th, the space 

W l ^\T K ) := {u\ Tk : u G W l ^ \T h )}, 

and in the same manner we define the spaces W l,p<y '\T z ) and W 1 '^' (T e ) , for any z G Mh and 
e ££ h . 

For each face e G 5^, e C Ti n t we denote by k + and k~ its neighboring elements. We write 
u + , v~ to denote the outward normal unit vectors to the boundaries dn^, respectively. The jump 
of a function u G W l,p ^'\Th) and the average of a vector-valued function <p £ (W 1 ' p ^'^(7fe))^ v , 
with traces , from fc 1 * 1 are, respectively, defined as the vectors 

[«] := u + ^ + + u~^~ and {0} := 



2 

Let h : <9Q U Ti n t ->la piecewise constant function define by 

h(x) = diam(e) if x G e, 

where e E £h- 

We consider the following seminorms on VF 1,p (')(7h), 

-i 

I«lwi*(-)f7fc) = II Vu IIlp(-)(Q) + IIM h7R llLP(-)(r lnt )> 

-1 

and the following local seminorm 

\ u \w^(-)(T K ) = ll Vu llLf(')(T K ) + Yl IIM h7W llxp(0(e). 

eCT K 

for any re G 7/j- Similarly, we define the seminorms \u\ w i, p (-)(t z ) an d \ u \w 1 'P<--'){T e ) f° r an y z ^ ■Mi 
and e G £h ■ 

Lemma 3.6. For aZZp: [1, oo) — > R, i/iere exisi a constant C, independent of h such that, 
\Du\(Q) < C\u\ w i, pi . ){Th) Vu G W 1A -\T h ), Vh G (0, 1]. 

Proof. For all u G W l ^\T h ), we have that 

|Du|(n)< I \\7u\dx+ I \M\ds. 

Thus, by Holder inequality, Proposition 2.1 (5) and the Hypothesis 3.2, there exists a constant 
C depending only of |fi|, p\ and p2 such that 

\Du\(tt) < C (\\Vu\\ LP , Kn) + ||h^) MlUo(rw)) • 

The proof is now complete. □ 

Lemma 3.7. Let (Th)he(o,i] be a family of partitions ofVt. Then, for each function p, q: — > [1, oo), 
there exists a constant C > independent of h, such that for any k G Th 

jv jy_ 

IMI LP (.) (K ) < Ch P K + q \\u\\ Lq , ){K) Vu G S k (%), Vh G (0, 1]. 
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Proof. Let k S 7~h, k its corresponding reference element and -F K : k — > k the associated affine 
mapping. We set J = \ det(DF K )\. Using Hypothesis 3.2, we have C~ 1 h K v < J < Ch^ , for some 
constant C which is independent of k. Let K > 0, then we have 



Thus, 

KChZrWMvuw < hoF K \\ LpaFKi . m . 

Using that /i K < 1, we obtain 

(3.H) \W\lv(-)(k) < (Ch%) 1/p+ \\uoF K \\ LpoFK( . Kk) . 

Similarly, we have 

(3-12) \\uoF K \\ LM , {k) < {Ch- N ) l li-\\u\\ Lq( , {K) . 

As on a finite dimensional space all the norms are equivalent, we have that there exists a 
constant C depending only on N and k such that, 

(3.13) \\uo F K \\ LpoFK( . ){k) < C\\uo F K \\ LP2(k) <C\\uoF K \\ Lqi{ f.) < C\\u o F\\ Lq( . ){k) , 

where in the first and last inequalities we are using Theorem 2.2. 

Finally, by (3.11)-(3.13) we arrive at the desired result. □ 

Lemma 3.8. If p is log-Holder continuous then, for any e £ D d£l and z £ Mh De we have 
that, 

i_ 

(3.14) IMlL p (->( e ) ^ Ch * P ~ \\u\\ LP (. ){Tz) Vn G S k (T h ), 
where C = C(pi,p 2 > N, 0, Q og ). 

Proof. Let k G Th such that e C k. Let F K and k be as in the proof of Lemma 3.7 and let 
e = F~\e). 
Then, 



k J I ^ V k 



Hence, 



1 II 

{C- l K)—) 



h N/p(x) 
"•re 



By using Theorem 2.2 and that all the norms are equivalent, we have 

\\u°Fk\\lp°Fk(-)(&) < C\\U0 F K \\ LP2 ^) < C\\U0 F K \\ L l(gy 

On the other hand, by the local inverse estimate in [6, page 837 ], we have 

\\ u ° F *\W(e) < C \\ u ° F K.\\v-(fiy 
By using again Theorem 2.2, we obtain 

\\ u ° ^relli^re) < C\\u o F K \\ LpoFK (-)( k y 

By using all the inequalities and the definition of the Luxembourg norm, we arrive at 
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1 U 



h N/p(x) 



< c 

LP(')( e ) 



h N/p(x) 

ll'K 



Lp()(k) 



N(p_-p + ) 
P-P+ 



Finally, we obtain 

i 

\\ h K p(x) u\\ L P(-)(e) ^ Ch K \\u\\lp(')( k ), 

By Remark 2.11, we get 

_L_ TV-Sr 

\\h K P^u\\ LP( . ){e) < Ce p i \\u\\ LP (. ){K y 

Now, inequality (3.14) follows immediately using Proposition 2.1(3) and the fact that h z ~ 
K. □ 

The next result establishes the existence of the local projector operator (for the proof see 
Subsection 3.1 of [6]). 

Lemma 3.9. For all z S Mh there exists a linear map tt z : BV(Cl) — > K such that 

\\u - tt z (u)\\ l i {Tz) < Ch z \Du\{T z ) Vn G BV(Sl) 
where C is a constant independent of h and z. 

4. The reconstruction operator Qh 

In many Discontinuous Galerkin problems one uses a priori bounds in order to prove the 
Poincare inequality for the discrete space. In order to prove these inequalities it is required to 
use a reconstruction operator. In this section we define, as in [6], a family of quasi-interpolant 
operators and prove some error estimates depending on the mesh size. These results are more 
general than the ones in [6], because we prove bounds in the variable p(x)— norm. On the other 
hand these results are less general than the ones in [4] in the sense that they only cover the case 
of the finite dimensional space S k (Th)- This last restriction comes from the fact that in Lemma 
3.7 we need to use the equivalence of the norms in the space of polynomials. 

In order to prove these error estimates we strongly use Proposition 2.11. This is the reason 
why we need p to be log-Holder continuous. 

Now, we define and study the reconstruction operator. For each h 6 (0, 1], let 

Q h : S*Cft)->-W 1,00 (fi) 

be the linear operator defined by 

Qh(u) = ^2 

z£Af h 

where A 2 is the standard P 1 nodal basis function associated with the vertex z on the mesh 7ft.. 

In the next theorem, we give some local estimates of the L q ^\n) and L q ('\e) norms in terms 
of the W 1 '^ (T K ) seminorm and h. 
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Theorem 4.1. Zei p, q: Q — > [l,oo) be log-Holder continuous in £1. Then, the operator Qh 
satisfies 

(4.15) \\u - Qh(u)\\ L g(.)( K) <Ch K P \u\ w i, p {-)(t k ) Vk e Th, 

N-l N | ^ 

(4.16) II^-Q/iWIIl^CO^) <Ch e q ~ \u\ w i, P (-)(T e ) Vee£ h ndtt, 

(4.17) I|VQ/>(u)|| Lp( . )(k) ^Cl-uj^i^.)^) Vk G T h , 

for all u € S k (7~h) where C is a constant independent of h. 

Proof. We proceed in three steps. 

Step 1. We first show inequality (4.15). 

Fix k € Th- For z £ A/"h fl k, by using Proposition 2.1 (6) and Lemma 3.7, we get 

f—N 

\\u-TT z (u)\\ Lq (.) iTz) <C 2^ Hn-vr^n)!^!^,). 

{k':k'CT z } 

By Remark 3.4 and Proposition 2.11, we get 

\\ u ~ ^z(u)\\ Lq (. )(T ^ < Chi \\U - 1T z (u)\\ L i {Tz y 

Thus, by Lemma 3.9, we have 

\\u- Tr z (u)\\ Lq (. ){Tz) <Cht ( ||Vu|| L i(7v) + [ \l u l\ ds ) • 

Then, by using again Lemma 3.7 and Remark 3.3 we have 

N 

eCT z 



JL +1 ( -2L r 

(4.18) \\u- Tr z (u)\\ Lq( . ){T;i) <Cht ( h z P ~ ||Vu|| LP (.) (T2) + K N ^2 J\ 



(4.19) 



To estimate the second term, we use Holder inequality, obtaining 

M\ds <2\\ Mh^\\^. Ke) \\h^\\^ { . He) 

<C\\ lujh e 7W || £ p(.)(e)/*e II 1 \\ L p'(-) ( e ) • 

Now, by Proposition 2.1 (5), we have that 



(iV-l)(i-^) 

||l|lip'(-)(e) - 

Then, we obtain 



N(l — M 

U^\ds<C\\Mhe P{X> \\r P (-)^h z 



l£P(-)(e)"2 

Therefore, summing over all e C T z and using (4.18), we arrive at 



aL-JL+i 

<?_ P- 



(4.20) ||it - vr 2 (n)|| L9 (. )(Tz ) < C/i z Mw/i,p(-)(t z )- 

Now, as in the proof Theorem 3.1 in [6] and using Proposition 2.1 (6), we have the inequality 
(4.15). 

Step 2. We now show the inequality (4.16). 
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Fix e G £h n dQ and let z G Mh H e. By the inequality (3.14), 



\\ u ~ Kz(u)\\Li<.-)(e) ^ Ch » W u ~ ^MILaOCT*)- 

Again, following the lines in [6] and using that p and q are log-Holder continuous in Q, we arrive 
at the inequality (4.16). 

Step 3. Finally, we will show the inequality (4.17). 

Fix k G 7~h- First, since (\ z )zeAf h r\K is a partition of the unity in k we have that for any x G k 
VQ h u(x) - Vu(x) = (k z (u)-u(x))V\ z (x)+ ^ Vu(x)X z (x) 

z£j\f h nn zGAfhClK 

\\VQhu\\ LP (.) {K) < \\M u )- u )v**\\ia>(-)(k)+ II VuA ^IIlp(')( K ) + \\Vu\\ LP (. )iK) . 

zeMu<~\K zeAfh<~\K 

Now, using Hypothesis (3.2), we have that there exists a constant C\ such that |VA 2 | < Ci/i" 1 
in k, and by (4.20) we get, using Remark 3.3, 

Lp('~) (k) W 1 ^') (T K ) 

<(C + 1)\ 

zeMhC\K 

The proof is now complete. □ 

Our next aim is to prove some global estimates. To this end we will need some definitions. 
Definition 4.2. Let p : O — > [1, oo). Given q : CI —> [1, oo) and q < p* in 0, we define 

f l ?( 3; ) 

7 = sup l^M :xG /' 

Observe that < 7 < 1 and 7 = if p(x) > N for all x G and 7 = 1 if p(x) < iV and 
q(x) = p*(x) for all 

Definition 4.3. Let p : fi — >■ [1, 00). Given q : — > [1, 00) and g < in 0, we define 

Ip*(x) J 

Observe that < /3 < 1 and /3 = if p(x) > iV for all x G $7 and /3 = 1 if p(x) < N and 
g(x) = p*(x) for all 

Lemma 4.4. Zei p,}:!!^ [1, 00) be log-Holder continuous in f2. Let u G S k (Th) satisfy 

( 4 - 21 ) \ u \w^-)(T h ) ^ L 

Then, we have that 

• If P < q <P* in fl, then 

(4.22) / |u-Q h (u)| 5(a!) dx <Ch N ^\ 



(4.23) / \VQ h (u)\P^dx<C, 

Jn 
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• If P < Q < P* in £1, then 

(4.24) f \u - Q h (u)\ q{x) dS < Ch^-W-ft, 

Jan 

where C = C(pi,p2,Q,Ci og , N) and 7 and (3 are given in Definitions 4. 2 and 4-3, respectively. 
Proof. First observe that, by (4.15), we have 



/ J5L_J2L+l V 



and by (4.21), we get 

1 



AT — hg- 



f \u-Q h {u)\ q ^dx < 1 Vk £ T h . 



Then, by Proposition 2.11, 

N — hg- 

rn, p - i„, 

lwA>;K-)(T re ) - K 1 W^H^) 



i« - Qh(«)r (a:) ^ < ch K p - \u\ q w ^, )!T , < ch K '\ u \ q - 



for any k £ 7~h and x G k. Therefore, 

(4.25) f \u-Q h (u)\^Ux<Ch N ^\u\ q wlM , {TK) V K GT ft . 

On the other hand, by Remark 3.3, the number of k C T k is uniformly bounded in h. Using 
this fact and Proposition 2.1 (6), we have that 

( 4 - 26 ) MW>py ^ c E (liv«fe., w + IIHh^|li;, )(Knrint) ) • 

On the other hand, if we suppose that || Vu^c-w^ > h^ q ~ , by Proposition 2.11 (2), we have 
that 

( 4 - 27 ) II^UW) < ^llV«||g(o w - 



Arguing as before, if \\lu}h p \\ lp() ( K r\T int ) — ' we ^ ave 

( 4 - 28 ) HMh^ir L ;, )(Knrmi) < c||Mh^ir L ;, )(Knrmt) . 

Now, we take 
and 

Observe that 



A={ K eT h :\\Vu\\ LP , ){K) >h^-}, 
B = { K eT h : HMh^lU^nr^) > h^~) . 



(4.29) 



E iiv«i&., w < ZX < c 

k€A c k^A c 

EllNh^||-, (Knrmt) <^ 
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On the other hand, by hypothesis (4.21), we have that ||Vn|| iP (.)^ K ^ < 1 and then, for all 



(4.30) 



iiv^iis C o W <iiv<:, )(K) </jv<Wdx, 

IIHh^HS ( . )((tnriBt) < / IMI^h 1 -^*. 

•/ /SHI 7'n.+ 



Since each /c appears only in finitely many sets T K > we have, by (4.26)— (4.30), 

E MW>(T.) ^ (e iiv«iiS(o w + E iiv«iiS(o w ) 

KeTh VkgA ksA c / 

+^(EiiHh^e ( . )(Knrmt) + E iiNh^n-, )(Knrint) 
- c ( E / 1 Vu i p(x) dx + E / i M p^h 1 ^^) + 1 ) 

=C ( / I Vul*) dx+ f | M pWh 1 ^^) ds + l). 
Thus, by (4.21) and (4.25) we get, 

f \u-Q h (u)\ q W dx=J2 [ \u-Q h {u)\ q{x) dx <Ch N ^\ 

KG lh 

Lastly, using the same argument, (4.16) and (4.17), we get 

\u - Q h {u)\ q{ - x) dS < c/i(tf-i)(i-0 and / \VQ h {u)\ p{ - x) dx < C, 



I an 

where C is independent of h. □ 
The following corollary follows immediately 

Corollary 4.5. Let p, q : Q [1, oo) &e log-Holder continuous in Q. Then, for all u G S k (Th), 
we have, 

• If p < q <p* inVL, then 



(4-31) ^ \u - Q h (u)\*W dx <Ch N ^ max {|n|^, p( . )(Th) , H^,^ ) } , 

(4-32) f iVQfcHI'W^^CrmMllu^^j.M^., }. 

• If P < q < P* in fl, then 
(4-33) |« - Q fc (t*)|«W d5 < Cfc<"-Da-0 maxj^l^^^, |< 1>K . )(7 - fe) } • 

Where C = C(p±,p2, ^, Ci og , N) and 7 and (3 are given in Definitions 4-2 and 4-3, respectively. 
Proof. It follows by Lemma 4.4, taking v = u\u\^ ltP ^^.. □ 
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Remark 4.6. Under the same hypothesis of the last corollary, if 1 < q < p* in Q, we have that, 
for all u G S k (%), 

\\u - Qh(u)\\ Lq( . ){n) < Ch N ^-^\u\ WUv( . ){Th) and \\VQ h (u)\\ LP( . ){n) < C\u\ wl , H . ){Th) , 
where C = C(pi,p 2 , ^, Ci og , N). 

5. The lifting operator 
We begin this section by defining, as in [6] (see also [1]), the lifting operator, i.e. 
Definition 5.1. Let I > and R h : W^iT) -> S l (%) N defined as, 

f {R h {u),<j>) dx = - ! (H, {0}) dS G S\%) N ■ 

This operator appears in the first term of the discretized functional Ih- As we can see from the 
definition, this operator represents the contribution of the jumps to the distributional gradient. 
This is the reason why it is crucial to add this term in order to have the consistency of the 
method. 

We point out that this lifting operator was first used in [2] in order to describe the contributions 
of the jumps across the interelements of the computed solution on the (computed) gradient of 
the solution in a mixed formulation of Navier-Stokes equations. It was also used in [5] where a 
solid mathematical background for the method introduced in [2] was proposed. 

Now, we give a bound of the L p (') (fJ)-norm of Rh(u) in terms of the jumps of u in T^f. 

When p is constant the proof follows from an inf-sup condition. Since in our case, we are 
dealing with the Luxemburg norm, we can't prove the boundedness directly from the definition. 
We can prove this inf-sup condition, but we can not use it to prove the result. Instead, we find 
a local characterization of Rh in order to prove a local bound and then, we prove the global 
bound. 

We give first the local estimate. 
Lemma 5.2. There exists a constant C\ such that, for any k G Th, we have 

\\Rh(u)\\ LPi , {K) < C||h-W(*) M | U)(KnIW V n G W^XTu) Vft G (0, 1]. 

Proof. We proceed in two steps. 

Step 1. We first want to prove that, 



(5.34) \R h {u)\ < / IMIdS G T h 

where C is independent of k and h. 

We begin by observing that, by Hypothesis 3.1, there exists m = m(k,N) G N such that for 
each k G Th, 



in 



Rh(u)\ K oF K = y~] anpi(x), 

i=l 

where {y?i} is the standard nodal base of (P l ) N in the reference element k := F^ 1 ^). 
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Using the definition of Rh we have that for each 1 < j < m, 

in 



ft lib ft ft 

/ R h {u)vj o F~\x) dx = Y,<H <Pi° FkH^W ° F-\x) dx = -Y J MiVj o F~\x)} dS. 

On the other hand, if we change variables and we use Hypothesis 3.2 and the fact that \<fi(x) | < 1, 
we get 

tpi o F~ 1 (x)ipj o F~ l (x) dx = h% [ ipi(x)ipj(x) \ dei ' D J K dx = d {j 

J k tl K . 



with dij ~ 1. 
Therefore, 

. m 

= j^w'^2( D ~ 1 b) i (p i (x)dx, 

K 8=1 

where D = (<iy) and 6j = — / [puK^j ° F f 7 1 (x)}d5. 
Thus, using that < 1, we arrive at (5.34). 

Step 2. Now, we show that there exists a constant C\ such that, for any k € Th, we have 

\\R h (u)\\ LP , HK) < C\\h-^'^M\U HK nr mt ) Vn € W^\%) Vh € (0, 1]. 
By inequality (4.19), we have 

N(l — U 

Mlds^c/ie p ||H^ PW || iP (.)( e) . 
Thus, by Hypothesis 3.2 and (5.34), we have that 



l_ 

Now, take T = ]>] egK ||H/i e p ' (a0 || LP (. )(e) . Then, 



R h (u) 



< C / hZ N ^ x)/p - dx < ChW-'+fr-) < C 



T 

where in the last inequality we are using Proposition 2.11. 

The result follows now by Remark 3.3. □ 

Lemma 5.3. Let p : Q — >: [1, oo) be log-Holder continuous in Q. Then, there exists a constant 
C such that, 

\\R h {u)\\ LP , ){n) <C\\h- l /v'^iul\\ LP , ){Vmt) VueW l ^(T h ) Me (0,1]. 

Proof. First, if we assume that ||h _1 / p ^Mll (r int ) — 1> we can P rove using Lemma 5.2 and 
proceeding as in Lemma 4.4 that, 



\R h {u)\ p W dx < C. 

Finally, taking v = u ( Wh" 1 ^' ^ lu}\\ LP (.) ^ r . nt ^\ , we obtain the desired result. □ 
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6. Convergence of the method 

In this section we first prove the broken Poincare Sobolev inequality which is crucial to get 
compactness. We also prove the coercivity of the functional and we finally arrive at the proof 
of Theorem 1.1. 

Theorem 6.1. Let p : Q, — s> [1, +oo) be log— Holder continuous in 0,. There exists a constant C 
such that, 

\W - (u)n\\ LP * ( . )(n) < C\u\ w i M . ){Th) Vu£S k (T h ) V/jG (0,1], 
where (u)q = -py J^udx. In particular, 

\\u\\lp*(-)(u) ^ ^(IHIl 1 ^) + \u\ w i, P (-)(j- h )) Vu£S k (T h ) V/iG (0,1]. 
Proof. We begin by observing that 

\\u - («)o|| L p*(-)(n) ^ W u ~ Q/*Mllzj>*(-)(n) + \\Qh(u) - (Qh(u))n\\LP*(-)(n) + C\\Q h (u) - u\\ L i {n] . 
Then, using the Remark 4.6 and Theorem 2.8, we have 

II" - («)oJ| LP *c)(n) ^ c \ u \w x >p(-){%) Vu G sk (Th) v/j 6 (°> !]■ 
The proof is complete. □ 

Theorem 6.2. For each h 6 (0, 1], let Uh £ W l,p<y '\Th)- If there exists a constant C independent 
of h such that Ih(uh) < C for all h £ (0, 1], then 



sup ( ||u/i||z,i(n) + Whlw^p^HTh) ) 



< oo. 



Moreover, 



sup I \u h - UD^h 1 -^ dS <oo. 
efo.il Jan 



he(o,i] Jon 

Proof. Since Ih( u h) — C then, ||h _1//p ' ( ' a: )[[u/ l ]|| i p(.)( r . iit ^ < C. And, by Lemma 5.3 and Proposi- 
tion 2.1, we have 

/ \R h {u h )\ p ^dx<C. 
Jn 

Using the third inequality in Proposition 2.5 and the above inequality, we get 

f \R h (u h ) + Vu h \ p(x) dx > 2 1 " P2 f \Vu h \ pi - x) dx-C. 
Jn Jn 

Therefore, 

I h {u h ) + C>2 l - p2 I \Vu h \ p ^dx + f \u h - u D \ p ^h l - p ^ dS + f Wuhlf^h 1 -^ dS. 

Thus, as Ih(uh) < C, we obtain that luhlw 1 ^ ) (%) an< ^ / \ Uh ~ ud\ p( " x ^ 1 ~ p ^ dS are uniformly 

J on 

bounded. 
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Finally, by Friedrichs inequality for BV, Lemma 3.6, Holder inequality, Proposition 2.1 and 
the fact that h < 1 we have, 

IWIl^Q) <C y u h\w 1 'P<--')(T h ) + J \ u h\dS 

<C (\u h \ w iM-) { r h ) + / \u D \dS + ||K -^)h~ 1/p ' (:E) || LP (.)(r I5 )l|h 1/p ' (x) l! LP '(.)(r I3 ) 
V J r D 

<C (\ u h\w^)(T h ) + \ u D\dS + \\(u h - u D )h" 1/p ' (:c:) || Lp( . )(rD) ^ . 

This completes the proof. □ 

Lemma 6.3. Let p, s and t be functions satisfying (HI), (H2) and (H3) respectively. Let 
Uh € S k (Th) be under the conditions of Theorem 6.2. Then, there exist a sequence hj — > and 
a function u € W 1,p ^(f2) such that 

(6.35) Ufy u weakly* in BV(SX) 

(6.36) Vu hj + R h {u hj ) Vu weakly in L p ^(0), 

(6.37) u^. — > u strongly in L S ('\Q), 

(6.38) uy lj — > u strongly in L^' (dU) . 

Proof. We first observe that, by Theorem 6.2, we have that 

sup (IM^i^) + \uh\ w x,v{-)(r h )) < 00 v/l e (°) !]• 
fte(o,i] 

Then, the proofs of (6.35) and (6.36) follows by applying Theorem 5.2 in [6] with value p\, and 
using that Rh(uh) is bounded in L p ('\{l), see Lemma 5.3. 

We now prove (6.37). By (6.35) and the compactness of the embedding BV(Q) C L 1 (fi), 
there exists a subsequence of nf,,. , still denoted by Uh 3 , such that 

Uhj — y u in L 1 (0). 

Since H^Hli + \ u hj\w x 'P^{Th) ls bounded, by Theorem 6.1, ||it/y ||xj>*c-)(a) is bounded, and by 
Theorem 2.9, u € W l ^\Vt) C L p *( )(fi). Therefore, using Theorem 2.3, we obtain that 

(6.39) u hj ^u in L s ^(fl), 

for all s satisfying (H2). 

Finally, we prove (6.38). We begin by observing that, by Corollary 4.5, 

(6.40) \\u h - Qh{uh)\\Lp(')(U) ~> as h -4 

and {Qh(uh)} is bounded in W 1,p ^'' (Q). Then, there exists v 6 W 1,p (')(fi) and subsequence 
{Qh,-(«fc 7 -)} such that Q hl {u hi ) -± v weakly in W lj, (-)(£i). Therefore, by (6.39) and (6.40), 



3 v 3 

V = u. 



Using Theorem 2.10, 

(6.41) Q h . (u hj ) -4 u strongly in L*« (dfi) . 

Now, taking t : Q — >■ [1, oo) log-Holder with t <i < p* and by Corollary 4.5, we get Q/^ (it/^ ) — 
Uh- — > strongly in L*W(9fl). Therefore, it^. — >• it in L*(')(9f2). □ 
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Before proving the convergence of the minimizers, we need an auxiliary lemma. It is in this 
step where we need more regularity of the boundary data. 

Lemma 6.4. Let h € (0,1], and p : Q, — > (l,oo) satisfying (HI). Assume that ud G W 2,P2 (Q) 
and let v € W 2 ' P2 (VL) n A then, there exists £ V (Th)> such that 

\\ v h - ->■ ash^O, 

and 

Ih{vh) -> I{v) as h ->• 0. 

Proof. Since p is log-Holder, we have that C°°(J)) are dense in W 1,p (')(Jl). Then the first part 
follows by standard approximation theory, see in Theorem 3.1.5 [9]. 
Moreover, Vh satisfies 

(6-42) \\v - v h \\ LP2 (dK) < C\\v - v h \\ w i, P2[K) < Ch K \\D 2 v\\ L v 2 ( K ) 



for each k 6 Th- Using Remark 3.4 and summing over all e £ dfl, we have 
(6.43) I \v- v h \ P2 h l ~ P2 ds < Ch\\D 2 vf 2 P2{n) . 

In addition, by Holder inequality and since h < 1, we have 

\v - v^h 1 -^ ds < C\\\v - v^h^-^/^W^,^. 



Since 



(\v - Vh \P( x )h ( - 1 - p2 ^ x ^ P2 ) P2 /P^ ds= [ \v- v h \ P2 h^- p2 ^ ds, 
an Jan 



then, by (6.43) 

\v - Vhl^h 1 -^ ds ^ as/t^O. 



I an 

Since v h £ W lj> U(Sl) then = and R h (v h ) = 0. Finally, using (6.42) and Theorem 2.6, we 
obtain the desired result. □ 

Now, we are in a condition to prove Theorem 1.1. 

Proof of Theorem 1.1. By Lemma 6.4 there exist Wh G C/ 1 (7/ l ) such that Wh — > urj strongly 
in Ty 1,p (')(fi) and Ih(wh) — > I{ud)- Therefore, since Ih(uh) < T\(wh), we have that Ih(uh) is 
bounded. 

Then, by Theorem 6.2 and Lemma 6.3 there exist a subsequence u^. and it 6 W 1,p ^ (O) such 
that 

it^ it weakly* in .BV(Jl), 
Vu^.+lJ^tt/yJ-^Vu weakly in ) (ft), 

Uhj — > it strongly in L s ^(ft), Vs satisfying (i?2) 
u/j. — » it strongly in L*^(9fi), Vt satisfying (H3). 
On the other hand, since the penalty term, 



/ h 1 "^ -u D | p dS 
•/iv, 
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is bounded, we have that 

\\u - Ud\\lp( )(T d ) - H n ~~ Uh j Wlp<--)(T d ) + W u hj - u d\\lp(-)(V d ) ~^ °- 

Then u € A. 

Taking s = q and t = r in (6.44), by Proposition 2.6 we have 



I{u) < liminf 

(6.45) 

< lim inf I h . (u hj ) < lim sup I h . (u h 



| (\Vu h] + R h (u hj )fW + [« hj - dx + J^ 



u h /W dS 



Now, we want to prove that u is the minimizer of /. Let v € Ar\W 2,P2 (Q), and let Vh € U l {J~h) 
as in Lemma 6.4. Then Ih{vh) ~^ I( v )- Therefore, by (6.45) 

(6.46) I(u) < lim inf I h (u h ) < lim sup I h . (u hj ) < lim hAv hj ) = I(v). 

j^oo j->oc J^°° 

Now, let w € A, then for any e > there exists v € .An VF 2 ' P2 (f2) such that \\v — Hlw^pOrn) < 
e. By Theorem 2.6 we have that I(v) < I(w) + e, therefore by (6.46) 

I(u) < lim inf I h . (u hj ) < lim sup I h . (u hj ) < I(w) + e. 

Taking e — > 0, we get 

J(u) < lim inf I h . (u hj ) < lim sup I h . (u hj ) < I(w) V w £ A. 

Therefore I(u) < I{w). 

Moreover, taking w = u, we have that all the inequalities in (6.45) are equalities, therefore 
we have that J/, 5 . (it/^. ) — > I(u). Then 

and using Lemma 5.3 we have that Rh{uhj) — > 0. This fact and (6.44) imply that Vu^- — 1 Vu 
weakly in U>^{p). 

Since u is the unique minimizer of /, the whole sequence Uh converges to u. 
Finally, since 



x Ux, 



Vu h + R h {u h ) Vu weakly in L P ^(Q) and f (\Vu h + R h (u h )\ p{x) dx -)■ f |Vu| p( 

by Proposition 2.6, Vii/, + Rh(uh) —> Vti strongly in L p ('\tl). Therefore, since Rh(uh) — > 
strongly in L p (')(f2), we get that Vuh —¥ Vtt strongly in L p (')(r2). □ 

7. The Continuous Galerkin Method 

In order to make a complete study of this problem, we prove the convergence of the Continuous 
Galerkin finite element method for our problem. In the next section, we make a comparison of 
the two methods in an example. 

For simplicity, we take the following functional: 

f ( |Vd p ^ \u-^\ , 
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with q(x) satisfying (H2). Then, since the functional I is strictly convex and coercive in A there 
exists a unique minimizer of the problem. 

We take now a partition of f2 as in Hypothesis 3.1 and the usual conforming subspace U k {Th) 
of This subspace consists of all continuous functions such that they are polynomials 

of degree at most k in each k G Th • We assume that for some h', ud G U k (Th>) (this assumption 
replaces the one on Lemma 6.4). 

Let now h < h! and 

= {v h G U k (T h ) :v h = u D on OQ}. 
For simplicity, we may assume that hi = 1. 

Remark 7.1. Let LT/j : Cq°(J7) — > be the interpolant mapping defined in Theorem 3.1.5 in [9]. 
Then, we have that 

\\Rh4> - 0llwri,P(-)(fi) - CH n /i^ ~ 011^^2(0) -> 0, 

for any ^> G Co°(fi). We also have, by the continuity of I, that I(ILh(f> + ud) — >• !{<P + Kfl) as 
/» -> 0. 

By the strict convexity of /, for each /i G (0, 1] there exists a function it^ G V£ such that 
is a minimizer in of /. 

Now we prove the main result of this section. 

Theorem 7.2. The sequence {u^} converges to u strongly in where u is the unique 

minimizer of I. 

Proof. Since {I(uh)} is uniformly bounded, there exists a subsequence of {uh} (still denoted by 
{uh}) such that 

(7.47) u h ->> u weakly in W 1,p(-) (n), 

(7.48) n fe -> u strongly in L p ^(n). 

As in the proof of Theorem 1.1, using Remark 7.1 instead of Lemma 6.4, we can prove that 
u is the minimizer of I and that I{uh) — > I(u) as h — > 0. By the convexity of 7 and (7.47), we 
have that 

(7.49) f \Vu h \ p{x) dx^ [ |Vn| p ^ dx as h -> 0. 

Then, by (7.48) and (7.49), using Proposition 2.6 (3), we have that {u^} converges to u strongly 
in W x ' p ^{n). 

Observe that, as in Theorem 1.1, we can conclude that the whole sequence {uh} converges to 
u strongly in W 1 '^^). □ 

8. One dimensional example 

In this section, we give an example in one dimension. Our idea is to compare the Continuous 
Galerkin Finite Element Method (CGFEM) versus the Discontinuous Galerkin Finite Element 
Method (DGFEM). We will see in the following example that, if the function p attains values 
close to one, our method converges faster to the solution. 
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Let Q = (-1, 1), < e, a < 1 and p : [-1, 1] ->• [1, 2] given by 



p(x) 



l-£ 

a 



scl + 1 + e if Ixl < a, 



2 if a < |x| < 1. 

For this function p(x) and for a given B > 0, we study the following problem, 



(8.50) 



'(K(x)p^)-V(x))' = in (-1,1), 
u(l) = --u(-l) = 5. 



We begin by observing that, since the operator is strictly monotone, we have a unique solution 

i 

of (8.50). Moreover, the solution satisfies |u'(x)| = CW*) -1 ) for some constant C > 0. Therefore, 
| it' (a;) | > max{C, C 1 ^} and, using that u € C l,a {[— 1, 1]), we have that it' does not change sign. 
Then, since u(l) > u(— 1) we obtain that u'(x) > 0. 
Thus, 

(8.51) u(x) = C(x + 1) - B if — 1 < a? < — o, 

(8.52) u(x) = C(a; - 1) + 5 if a < x < 1. 



Since p is even, we have that u is odd, so u(0) = 0. Therefore, u(x) = Jq C^- 1 ds for all 
x € [—a, a]. 

On the other hand, since the derivative of u at zero has modulus C 1 ^ 6 , if C > 1 we have 

lim \u'(0)\ = +oo. 

£->0 

This is reasonable since we expect to have big derivative when p approaches the value one. 



From now on, we take e = a = .01 and since it is easier to get B from C, we impose C = 1.3. 
Then B = ^ 13^^ ds ~ 1.03 10 6 . Observe that in this case, \u'{0)\ = 1.3 100 ~ 2.4 x 10 11 . 

Now, we find the corresponding solution for the CGFEM and the DGFEM. In both cases 
we take a uniform partition of [—1, 1] in n subintervals with size 2/n and k = 1. We use the 
trapezoidal numerical quadrature to compute the integrals appearing in the discrete functionals. 
The analysis of these integration errors falls beyond the scope of this paper. 

Observe that for the continuous method, we impose the boundary conditions and then, the 
space where we find minimizers has dimension n — 2. For the Discontinuous method, since 
we do not impose conditions on the boundary, and the number of nodal basis are 2n — 2, we 
are minimizing in a space of this dimension. Therefore, to make a comparison between both 
methods, we compare the discrete problem for the DGFEM in n— intervals with the CGFEM in 
2n— intervals. 

We want to mention that, in order to find minimizers of both discrete problems, we use a 
BFGS Quasi-Newton method (see [19, 25]). 

In the next two figures, we first plot the solution versus the approximation using the DGFEM 
and CGFEM for the case n = 41. The second figure is the graphic of the function p(x). 
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-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 




-1 -0.8 -0.6 -0.4 -0.0 0.2 0.4 0.6 0. 



Note that, when we use the CGFEM the discrete solution is close to the function y = x which 
is a solution of (8.50) with p = 2, that means that this method needs a smaller step in order to 
see the points where p is close to one. 

In the following figure we can see that, the minimizers of the continuous methods are far from 
the solution even for n = 150 (300-intervals). We need n = 200 (400-intervals) to arrive at a 
good approximation of u. 
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